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Preface 



What follows are my lecture notes for a first course in differential equations, 
taught at the Hong Kong University of Science and Technology. Included in 
these notes are links to short tutorial videos posted on YouTube. 

Much of the material of Chapters 2-6 and 8 has been adapted from the 
widely used textbook “Elementary differential equations and boundary value 
problems” by Boyce & DiPrima (John Wiley & Sons, Inc., Seventh Edition, 
©2001). Many of the examples presented in these notes may be found in this 
book. The material of Chapter 7 is adapted from the textbook “Nonlinear 
dynamics and chaos” by Steven H. Strogatz (Perseus Publishing, ©1994). 

All web surfers are welcome to download these notes, watch the YouTube 
videos, and to use the notes and videos freely for teaching and learning. An 
associated free review book with links to YouTube videos is also available from 
the ebook publisher bookboon.com. I welcome any comments, suggestions or 
corrections sent by email to jeffrey.chasnov@ust.hk. Links to my website, these 
lecture notes, my YouTube page, and the free ebook from bookboon.com are 
given below. 

Homepage: 

http : //www .math.ust . hk/~machas 



YouTube: 

https : //www. youtube . com/ user / j chasnov 



Lecture notes: 

http : //www. math. ust . hk/~machas/dif f erential- equations .pdf 
Bookboon: 

http : / /bookboon . com/ en/diff erential- equations -with- youtube- examples- ebook 
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Chapter 0 

A short mathematical 
review 



A basic understanding of calculus is required to undertake a study of differential 
equations. This zero chapter presents a short review. 



0.1 The trigonometric functions 



The Pythagorean trigonometric identity is 

sin x + cos x = 1, 
and the addition theorems are 

sin(x + y) = sin(a;) cos (y) + cos(x) sin(y), 
cos(a: + y) = cos(x) cos (y) — sin(a;) sin(y). 

Also, the values of sin a; in the first quadrant can be remembered by the rule of 
quarters, with 0° = 0, 30° = 7r/6, 45° = 7t/4, 60° = 7t/3, 90° = n/2: 




The following symmetry properties are also useful: 

sin(7r/2 — x) = cos a:, cos(7r/2 — x) = sin a:; 



and 



sin(— x) = — sin(x), cos(— x) = cos(a;). 



0.2 The exponential function and the natural 
logarithm 

The transcendental number e, approximately 2.71828, is defined as 

e = lim 

n—>o o 




l 
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The exponential function exp (x) = e x and natural logarithm In x are inverse 
functions satisfying 

e lnx =x, lne x = x. 

The usual rules of exponents apply: 

e x e y = e x+y ^ ^^y = & x-y , ^xy = e , P x 

The corresponding rules for the logarithmic function are 

In (xy) = In x + In y, In (x/y) = In x — In y, In x p = p In x. 



0.3 Definition of the derivative 



The derivative of the function y = /( x), denoted as f'(x) or dy/dx, is defined 
as the slope of the tangent line to the curve y = f(x) at the point (x,y). This 
slope is obtained by a limit, and is defined as 



/'(*) 



lim + h )~ f( x ) 

h — >-0 h 



(1) 



0.4 Differentiating a combination of functions 

0.4.1 The sum or difference rule 

The derivative of the sum of f{x) and g(x) is 

C f + 9)' = f' + 9 '■ 

Similarly, the derivative of the difference is 

(. f-9)' = f'-9 '■ 



0.4.2 The product rule 

The derivative of the product of f(x) and g{x) is 

( fa Y = f'g + fg', 

and should be memorized as “the derivative of the first times the second plus 
the first times the derivative of the second.” 



0.4.3 The quotient rule 

The derivative of the quotient of f{x) and g(x) is 

(f\ _ f'g- f g' 
\g) g 2 



and should be memorized as “the derivative of the top times the bottom minus 
the top times the derivative of the bottom over the bottom squared.” 
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0.4.4 The chain rule 

The derivative of the composition of f(x) and g(x) is 

(/(S'O))) = f'( 9 (x)) ■ g'{x), 

and should be memorized as “the derivative of the outside times the derivative 
of the inside.” 



0.5 Differentiating elementary functions 

0.5.1 The power rule 

The derivative of a power of x is given by 

-^-x p = px p ~ l . 
ax 

0.5.2 Trigonometric functions 

The derivatives of sin x and cos x are 

(sin a;/ = cos a;, (cos a:) 7 = — sin a;. 

We thus say that “the derivative of sine is cosine,” and “the derivative of cosine 
is minus sine.” Notice that the second derivatives satisfy 

(sin x)" = — sin x, (cos x)" = — cos x. 



0.5.3 Exponential and natural logarithm functions 

The derivative of e x and In a; are 

{e x )' = e x , (ln®V = -. 

x 

0.6 Definition of the integral 

The definite integral of a function f(x) > 0 from x = a to b (b > a) is defined 
as the area bounded by the vertical lines x = a, x = b, the x-axis and the curve 
y = f(x). This “area under the curve” is obtained by a limit. First, the area is 
approximated by a sum of rectangle areas. Second, the integral is defined to be 
the limit of the rectangle areas as the width of each individual rectangle goes to 
zero and the number of rectangles goes to infinity. This resulting infinite sum 
is called a Riemann Sum , and we define 

,6 N 

/ f(x)dx = lim X f(a + (n-l)h)-h, (2) 

71=1 

where N = (b — a)/h is the number of terms in the sum. The symbols on the 
left-hand-side of ([2]) are read as “the integral from a to b of f of x dee x.” The 
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Riemann Sum definition is extended to all values of a and b and for all values 
of f(x) (positive and negative). Accordingly, 

pa pb pb pb 

/ f(x)dx = — f{x)dx and / (—f(x))dx = — / f(x)dx. 

J b J a J a J a 

Also, if a < b < c, then 

[ f{x)dx = f f(x)dx + f f(x)dx, 

J a J a J b 

which states (when f(x) > 0) that the total area equals the sum of its parts. 



0.7 The fundamental theorem of calculus 



view tutorial 

Using the definition of the derivative, we differentiate the following integral: 



d 

dx 



f(s)ds 



lim ,C h md» fo)d» 

h — >-0 h 



lim 

h->- o 



lim 

h^O 



i: +h f(s)ds 

h 

hf (x) 
h 



f{x). 



This result is called the fundamental theorem of calculus, and provides a con- 
nection between differentiation and integration. 

The fundamental theorem teaches us how to integrate functions. Let F{x) 
be a function such that F'(x ) = f{x). We say that F(x) is an antiderivative of 
f(x). Then from the fundamental theorem and the fact that the derivative of a 
constant equals zero, 

F(x) = f f(s)ds + c. 



Now, F(a) = c and F(b) = f(s)ds + F(a). Therefore, the fundamental 
theorem shows us how to integrate a function f{x) provided we can find its 
antiderivative: 



f f(s)ds=F(b)-F(a). 

J a 



(3) 



Unfortunately, finding antiderivatives is much harder than finding derivatives, 
and indeed, most complicated functions cannot be integrated analytically. 

We can also derive the very important result J§ directly from the definition 
of the derivative ([l]) and the definite integral ([2]). We will see it is convenient 
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to choose the same h in both limits. With F'(x) = f(x), we have 

[ f{s)ds = f F'(s)ds 
J a J a 

N 

= lim F' (a + (n — 1 )h) • h 



F{a + nh) — F(a + (n — 1 )h) 

= lim > H • h 

h^O ' h 



n—1 

N 



= lim X F{a + nh) — F(a+ (n — 1 )h). 



The last expression has an interesting structure. All the values of F(x) eval- 
uated at the points lying between the endpoints a and b cancel each other in 
consecutive terms. Only the value —F(a) survives when n = 1, and the value 
+F(b) when n = N, yielding again 



0.8 Definite and indefinite integrals 

The Riemann sum definition of an integral is called a definite integral. It is 
convenient to also define an indefinite integral by 

J f{x)dx = F(x), 

where F(x) is the antiderivative of f(x). 



0.9 Indefinite integrals of elementary functions 



From our known derivatives of elementary functions, we can determine some 
simple indefinite integrals. The power rule gives us 



/ 



x n dx 



x n+l 

n + 1 



+ c ) 



n — 1. 



When n = — 1, and x is positive, we have 



— dx = In x + c. 
x 



If x is negative, using the chain rule we have 



_d_ 

dx 



In (—x) 



1 

x 



Therefore, since 

, , f —x if x < 0; 

|a; l = \* if ie > 0, 

we can generalize our indefinite integral to strictly positive or strictly negative 
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Trigonometric functions can also be integrated: 

J cos xdx = sin x + c, J sin xdx = — cos x + c. 
Easily proved identities are an addition rule: 

J (f( x ) + g{x))dx = J f{x)dx + J g(x)dx; 

and multiplication by a constant: 

J Af(x)dx = A J f{x)dx. 

This permits integration of functions such as 

f . o _ X 3 7x 2 

/ ( x 2 + 7x + 2)dx = — H — - — (-2 x + c, 

J o 2 



and 



J (5 cos x + sin x)dx = 5 sin x — cos x + c. 



0.10 Substitution 

More complicated functions can be integrated using the chain rule. Since 

^/(sW) = f'(g{x )) • g'{x), 

we have 

J f(g(x)) ■ g\x)dx = f(g(x)) + c. 

This integration formula is usually implemented by letting y = g{x). Then one 
writes dy = g'{x)dx to obtain 

J f'{g{x))g'{x)dx = J f'(y)dy 

= f{y) + c 
= f(g(x)) +c. 



0.11 Integration by parts 

Another integration technique makes use of the product rule for differentiation. 
Since 

ifg)' = f'g + fg '> 

we have 

f'g = {fg)' - fg'- 

Therefore, 

J f'{x)g{x)dx = f(x)g(x) - j f{x)g'{x)dx. 
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Commonly, the above integral is done by writing 



u = g(x) 
du = g' (x)dx 

Then, the formula to be memorized is 



dv = f'(x)dx 
v = f(x). 



I 



udv = uv — 



I 



vdu. 



0.12 Taylor series 



A Taylor series of a function f(x) about a point x = a is a power series rep- 
resentation of f(x) developed so that all the derivatives of f(x) at a match all 
the derivatives of the power series. Without worrying about convergence here, 
we have 

f\ x ) = /( a ) + / {a){x — o) H — [x — a) + , (x - a) +.... 

Notice that the first term in the power series matches /(a), all other terms 
vanishing, the second term matches all other terms vanishing, etc. Com- 

monly, the Taylor series is developed with a = 0. We will also make use of the 
Taylor series in a slightly different form, with x = x* + e and a = x* : 

f{x* + e) = f{x t ) + /'( x*)e + f ^*h 2 + - i^l e 3 + .... 

Another way to view this series is that of g(e) = f(x* + e), expanded about 
e = 0. 

Taylor series that are commonly used include 



= 1 



sin x = x — 



X X 

:+ 2! + ¥ 

y.3 ry, 5 



3! 



X 

5 !" 

A 



. X X 
CO S X = l-- + ~- 



1 



1 + X 



= 1 — X 



X 2 - 



r 2 3 

In (1 + x) = x - — + — - 



for \x\ < 1, 

, for lad < 1. 



A Taylor series of a function of several variables can also be developed. Here, 
all partial derivatives of f(x, y ) at (a, b ) match all the partial derivatives of the 
power series. With the notation 



fx = 



Of 

dx ’ 



we have 



fy ~ 



0J_ 

Oy ’ 



a 2 / 

dx 2 ’ 



0 2 f 

dxdy’ 



dy 2 ’ 



etc., 



f(x,y) = f(a,b) + f x (a, b)(x - a) + f v (a,b)(y - b) 

+ ( fxx(a,b){x - a) 2 +2f xy (a, b) (x - a) (y - b) + f yy {a,b)(y - b) 2 ) + ... 
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0.13 Complex numbers 

view tutorial: Complex Numbers 
view tutorial: Complex Exponential Function 
We define the imaginary number i to be one of the two numbers that satisfies 
the rule (i) 2 = — 1, the other number being —i. Formally, we write i = \j — 1. A 
complex number 2 is written as 



z = x + iy, 

where x and y are real numbers. We call x the real part of z and y the imaginary 
part and write 

x = Re z, y = Im 

Two complex numbers are equal if and only if their real and imaginary parts 
are equal. 

The complex conjugate of z = x + iy, denoted as z , is defined as 

z = x — iy. 



Using z and z, we have 



Re 2 = — (z + z) , Im z = — (z — z) . 
2 2 1 



Furthermore, 



zz = ( x + iy)(x — iy) 

2 -2 2 
= x — i y 

2 , 2 
= X +y ; 

and we define the absolute value of z, also called the modulus of z, by 

1*1 = (* 2) 1/2 
= y/x 2 + y 2 . 

We can add, subtract, multiply and divide complex numbers to get new 
complex numbers. With z = x + iy and w = s + it, and x, y, s, t real numbers, 
we have 



z + w = (x + s) + i(y + t); z — w = (x — s) + i(y — t); 

zw = [x + iy)(s + it) 

= (, xs — yt) + i(xt + ys)\ 

z zw 
w ww 

(x + iy){s — it) 
s 2 + 1 2 

_ ( xs + yt ) .( ys-xt ) 

s 2 +t 2 +l s 2 +f 2 ■ 
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Furthermore, 

\zw\ = \J ( xs — yt ) 2 + ( xt + ys ) 2 
= / (a; 2 + y 2 )(s 2 + t 2 ) 

= |z|M; 



and 



zw = (xs — yt) — i(xt + ys) 
= (x- iy) (s - it) 

= zuu 



Similarly 

1 - M TE\ - 1 

w \w\ w w 

Also, z + w = z + w. However, \z + w\ < \z\ + |ie|, a theorem known as the 
triangle inequality. 

It is especially interesting and useful to consider the exponential function of 
an imaginary argument. Using the Taylor series expansion of an exponential 
function, we have 






= 1 + 



= 1 - 



m - 

e 2 



(iO) 2 (iO) 3 (id) 4 (i6») E 



2 ! 

2! ' 4! ' 
= cos 6 + i sin 0. 



3! 



4! 

+ *( 9- 



03 

3! 



5! 

9 5 

+ 5! + " 



Therefore, we have 

cos 9 = Re e l ° , sin 9 = Irn e l 6 . 

Since cos7r = — 1 and sin7r = 0, we derive the celebrated Euler’s identity 



e in + 1 = 0, 



that links five fundamental numbers, 0, 1, i, e and 7r, using three basic mathe- 
matical operations, addition, multiplication and exponentiation, only once. 

Using the even property cos (—9) = cos 9 and the odd property sin (—9) = 
— sin 0, we also have 

e~‘ 6 = cos 0 — *sin0; 

and the identities for e™ and e~ ,e results in the frequently used expressions, 



cos 0 = 




sin0 



e id _ e -iB 



2 i 



The complex number z can be represented in the complex plane with Re 2 
as the x-axis and Im 2 as the y-axis. This leads to the polar representation of 
2 = x + iy: 

z = re ie , 

where r = | 2 | and tan0 = y/x. We define arg 2 = 0. Note that 0 is not unique, 
though it is conventional to choose the value such that — 7 r < 0 < 7r, and 0 = 0 
when r = 0. 
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Useful trigonometric relations can be derived using e ,e and properties of the 
exponential function. The addition law can be derived from 

e i{x+y) _ gix e iy . 

We have 



cos(a; + y) + i sin(a; + y) = (cos x + i sin x ) (cos y + i sin y) 

= (cos x cos y — sin x sin y) + i (sin x cos y + cos x sin y) ; 

yielding 

cos(a: + y) = cos x cos y — sin x sin y, sin(x + y) = sin x cos y + cos x sin y . 

De Moivre’s Theorem derives from e mS = {e l6 ) n , yielding the identity 

cos(nO) + isin(nd) = (cos0 + *sin0) n . 

For example, if n = 2, we derive 

cos 2 9 + i sin 29 = (cos 6 + i sin 9) 2 

= (cos 2 9 — sin 2 9) + 2 i cos 9 sin 9. 



Therefore, 

cos 29 = cos 2 9 — sin 2 9 , sin 29 = 2 cos 9 sin 9. 

With a little more manipulation using cos 2 9 + sin 2 9 = 1, we can derive 

o _ 1 + cos 29 . o „ 1 — cos 29 

cos 2 9= , sin -9= , 

which are useful formulas for determining 

J cos 2 9 d9 = ^ (2 9 + sin 29) + c, J sin 2 9 d9 = ^ (29 — sin 29) + c, 

from which follows 



o 



sin 2 9d9 = 



o 



cos 2 9 d9 = 7 r. 




Chapter 1 

Introduction to odes 



A differential equation is an equation for a function that relates the values of 
the function to the values of its derivatives. An ordinary differential equation 
(ode) is a differential equation for a function of a single variable, e.g., x(t), while 
a partial differential equation (pde) is a differential equation for a function of 
several variables, e.g., v{x,y,z,t). An ode contains ordinary derivatives and a 
pde contains partial derivatives. Typically, pde’s are much harder to solve than 
ode’s. 



1.1 The simplest type of differential equation 

view tutorial 

The simplest ordinary differential equations can be integrated directly by finding 
antiderivatives. These simplest odes have the form 



d n x 

~dtP 



G(t ), 



where the derivative of x = x[t) can be of any order, and the right-hand-side 
may depend only on the independent variable t. As an example, consider a mass 
falling under the influence of constant gravity, such as approximately found on 
the Earth’s surface. Newton’s law, F = ma, results in the equation 



d 2 x 

= -mg, 

where x is the height of the object above the ground, m is the mass of the 
object, and g = 9.8 meter/sec 2 is the constant gravitational acceleration. As 
Galileo suggested, the mass cancels from the equation, and 

d 2 x 



Here, the right-hand-side of the ode is a constant. The first integration, obtained 
by antidifferentiation, yields 



dx 

dt 



A-gt, 
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with A the first constant of integration; and the second integration yields 

x = B + At — ^ gt 2 , 

with B the second constant of integration. The two constants of integration A 
and B can then be determined from the initial conditions. If we know that the 
initial height of the mass is Xq, and the initial velocity is Vq , then the initial 
conditions are 

a;(0) = x 0 , < ^-(0)=v o . 

Substitution of these initial conditions into the equations for dx/dt and x allows 
us to solve for A and B. The unique solution that satisfies both the ode and 
the initial conditions is given by 

x{t) = x 0 + v 0 t - ^gt 2 . (1.1) 



For example, suppose we drop a ball off the top of a 50 meter building. How 
long will it take the ball to hit the ground? This question requires solution of 
(1.1) for the time T it takes for x(T) = 0, given x$ = 50 meter and vq = 0. 
Solving for T, 



T = 




sec 



3.2sec. 



Chapter 2 



First-order differential 
equations 

Reference: Boyce and DiPrima, Chapter 2 



The general first-order differential equation for the function y 
as 



dy_ 

dx 






y(x) is written 

(2.1) 



where /( x, y) can be any function of the independent variable x and the depen- 
dent variable y. We first show how to determine a numerical solution of this 
equation, and then learn techniques for solving analytically some special forms 
of (2.1), namely, separable and linear first-order equations. 



2.1 The Euler method 



view tutorial 



Although it is not always possible to find an analytical solution of (2.1) for 
y = y{x ), it is always possible to determine a unique numerical solution given 
an initial value y( xq) = yo, and provided f(x,y) is a well-behaved function. 
The differential equation (£3 gives us the slope f{xo,yo ) of the tangent line 
to the solution curve y = y(x) at the point ( Xo,yo ). With a small step size 
Ax, the initial condition ( Xo,yo ) can be marched forward in the x-coordinate 
to x = Xq + Ax, and along the tangent line using Euler’s method to obtain the 
y-coordinate 

y(x o + Ax) = y(x 0 ) + Ax/(x 0 , y 0 ). 



This solution (xo + Ax, yo + Ay) then becomes the new initial condition and is 
marched forward in the x-coordinate another Ax, and along the newly deter- 
mined tangent line. For small enough Ax, the numerical solution converges to 
the exact solution. 
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2.2 Separable equations 

view tutorial 

A first-order ode is separable if it can be written in the form 

= f(x), y(xo) = Vo, (2.2) 

where the function g(y) is independent of x and f(x) is independent of y. Inte- 
gration from xq to x results in 

fX nX 

/ g(y{x))y' {x)dx = / f(x)dx. 

j X 0 J Xq 

The integral on the left can be transformed by substituting u = y(x ), du = 
y' [x)dx, and changing the lower and upper limits of integration to y(x o) = yo 
and y(x) = y. Therefore, 



r 

'yo 



g[u)du = / f(x)dx, 

J Xn 



and since u is a dummy variable of integration, we can write this in the equivalent 
form 

[ g{y)dy= j f{x)dx. (2.3) 

Jyo dx 0 



A simpler procedure that also yields ( 2 . 3 ) is to treat dy/dx in (2.2) like a fraction. 
Multiplying (2.2) by dx results in 

g(y)dy = f(x)dx, 

which is a separated equation with all the dependent variables on the left-side, 
and all the independent variables on the right-side. Equation (2.3) then results 
directly upon integration. 



Example: Solve ^ 



dx ' \y= b with 1/(0) = 2. 

We first manipulate the differential equation to the form 



dy bo n 



(2*4) 



and then treat dy/dx as if it was a fraction to separate variables: 



dy 

3 -y 



-dx. 

2 



We integrate the right-side from the initial condition x = 0 to x and the left-side 
from the initial condition y(0) = 2 to y. Accordingly, 




(2.5) 
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x 



Figure 2.1: Solution of the following ode: ^ = |. 



The integrals in (2.5 1 need to be done. Note that y( x) < 3 for finite x or the 
integral on the left-side diverges. Therefore, 3 — y > 0 and integration yields 



— In (3 — y)] g = ^x] x Q , 

1h (3 — y) = -^x, 

3 -y = e -i*, 
y = 3-e~ h - x . 



Since this is our first nontrivial analytical solution, it is prudent to check our 
result. We do this by differentiating our solution: 

dy = - e ~h* 

dx 2 

= 2 ( 3 -^); 

and checking the initial conditions, y( 0) = 3 — e° = 2. Therefore, our solution 
satisfies both the original ode and the initial condition. 



Example: Solve ^ + \y = §, with y(0) = 4. 

This is the identical differential equation as before, but with different initial 
conditions. We will jump directly to the integration step: 



dy 

3 -y 



1 



dx. 



16 



CHAPTER 2. FIRST-ORDER ODES 



Now y{x) > 3, so that y — 3 > 0 and integration yields 



— In (j/ — 3)] ^ = *x] x 0 , 
ln(y-3) = -^x, 

y- 3 = e -i*, 

y = 3 + e“H 



The solution curves for a range of initial conditions is presented in Fig. 2.1 



All solutions have a horizontal asymptote at y = 3 at which dy/dx = 0. For 
2/(0) = 2 / 0 , the general solution can be shown to be y(x) = 3+(yo~3) exp(— a;/2). 



Example: Solve ^ = 2 g ° s 2 , with y(0) = —1. (i) For what values of 
x > 0 does the solution exist? (ii) For what value of x > 0 is y( x) 
maximum? 

Notice that the solution of the ode may not exist when y = —3/2, since dy/dx — » 
oo. We separate variables and integrate from initial conditions: 

(3 + 2 y)dy = 2 cos 2x dx 

/ y rx 

(3 + 2y)dy = 2 / cos 2a; da; 

-l Jo 

3// + ?/ 2 ] V _ 1 = sin 2x] “ 
y 2 + 3y + 2 — sin 2x = 0 

V± = ^h 3± vT + 4 sin 2cc] . 

Solving the quadratic equation for y has introduced a spurious solution that 
does not satisfy the initial conditions. We test: 

2 /± ( 0 ) = ^ [— 3 ± 1 ] = | 

Only the + root satisfies the initial condition, so that the unique solution to the 
ode and initial condition is 

y = ^ [— 3 + \J\ + 4 sin 2x] . (2.6) 

To determine (i) the values of x > 0 for which the solution exists, we require 

1 + 4sin2x > 0, 



or 

sin2a; > — (2.7) 

- 4 

Notice that at x = 0, we have sin2a; = 0; at x = 7t/4, we have sin2x = 1; 
at x = 7t/2, we have sin 2a: = 0; and at x = 37r/4, we have sin 2a: = —1 We 
therefore need to determine the value of x such that sin 2a: = —1/4, with x in 
the range 7r/2 < x < 37r/4. The solution to the ode will then exist for all x 
between zero and this value. 
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(3+2y) dy/dx = 2 cos 2x, y(0) = -1 




Figure 2.2: Solution of the following ode: (3 + 2 y)y' = 2cos2x, y(fi) = — 1. 



To solve sin2x = — 1/4 for x in the interval tt/ 2 < x < 37t/ 4, one needs to 



recall the definition of arcsin, or sin 
The inverse of the function 

/(x) = sinx, 



as found on a typical scientific calculator. 
— 7 t/2 < x < 7t/2 



is denoted by arcsin. The first solution with x > 0 of the equation sin 2x = — 1/4 
places 2x in the interval (7t,37t/2), so to invert this equation using the arcsine 
we need to apply the identity sin (n — x) = sinx, and rewrite sin2x = —1/4 as 
sin (7T — 2x) = —1/4. The solution of this equation may then be found by taking 
the arcsine, and is 

7T — 2x = arcsin (—1 /4) , 



or 

If .1 

x = — 7T + arcsm - 
2 \ 4 

Therefore the solution exists for 0 < x < (n + arcsin (1/4)) /2 = 1.6971..., 
where we have used a calculator value (computing in radians) to find arcsin(0.25) = 

0.2527 At the value (x, y) = (1.6971 . . . , —3/2), the solution curve ends and 

dy/dx becomes infinite. 

To determine (ii) the value of x at which y = y(x) is maximum, we examine 



(2.6 1 directly. The value of y will be maximum when sin2x takes its maximum 



value over the interval where the solution exists. This will be when 2x = 7r/2, 
or x = 7r/4 = 0.7854 



The graph of y = y{x) is shown in Fig. 2.2 



2.3 Linear equations 

view tutorial 
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The first-order linear differential equation (linear in y and its derivative) can be 
written in the form 



dy . . . . 

— +p{x)y = g(x ), 
dx 



( 2 . 8 ) 



with the initial condition y(x o) = yo . Linear first-order equations can be inte- 
nd). V 

-p(x)y 



grated using an integrating factor y(x). We multiply (2.8) by y(x), 

’ dy 



dx 



p{x) 

and try to determine y{x) so that 
y{x) 



%+vti)v 



= M(a ')g(x), 






Equation (2.9) then becomes 



— [y,(x)y] = y(x)g(x). 



(2.9) 



(2.10) 



( 2 . 11 ) 



Equation (2.11) is easily integrated using y(xo) = yo and y(x o) = y o'- 

g(x)y - y 0 y 0 = / y{x)g{x)dx , 

Jx 0 

1 



or 



V = 



H(x) 



Mo yo 



y(x)g(x)dx 



(2.12) 



It remains to determine y(x) from (2.10). Differentiating and expanding (2.10) 
yields 

dy dy, dy 



d~r +ppy = + 

dx dx dx 



and upon simplifying, 



dy 

dx 



= py. 



(2.13) 



Equation (2.13) is separable and can be integrated: 

r m 



r*=r*)*. 

1(10 ” J Xo 
C x 

In — = p(x)dx , 
Mo J xq 

y{x) = Mo ex P ( / p{x)dx 

\J X o 



Notice that since yo cancels out of (2.12), it is customary to assign yo = 1. The 
solution to (2.8) satisfying the initial condition y{x o) = yo is then commonly 
written as 



V = (do + J y(x)g{x)dx^j , 

y(x) = exp I / p(x)dz 

\J In 



with 



the integrating factor. This important result finds frequent use in applied math- 
ematics. 
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Example: Solve ^ + 2 y = e x , with y(0 ) = 3/4. 

Note that this equation is not separable. With p(x) = 2 and g{x) = e~ x , we 
have 



and 



p(x) = exp 

X2x 





— 2x 



— 2x 



— 2x 




+ 



e x dx 






1 

4 




Example: Solve gf — 2xy = x, with y(0) = 0. 

This equation is separable, and we solve it in two ways. First, using an inte- 
grating factor with p(x) = — 2x and g(x) = x: 



y(x) = exp ( —2 




and 



y = e 



xe x dx. 



The integral can be done by substitution with u — x 2 , du = 2xdx: 



rx o i rx 2 

/ xe~ x dx = - e~ u du 

Jo 2 J 0 




1 

2 



( 



)■ 




1 — e 




) 



Therefore, 
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Second, we integrate by separating variables: 



d V 0 

2 xy = x , 

ax 

| = x( 1 + 2 !) ), 



dy 

1 + 2 y 



= / xdx. 



1 



^ In (1 + 2 y) = + 

1 + 2 y = e x \ 

1 

V = 



(»'■ - 0 ■ 



The results from the two different solution methods are the same, and the choice 
of method is a personal preference. 



2.4 Applications 



2.4.1 Compound interest with deposits or withdrawals 

view tutorial 

The equation for the growth of an investment with continuous compounding 
of interest is a first-order differential equation. Let S(t) be the value of the 
investment at time t , and let r be the annual interest rate compounded after 
every time interval At. We can also include deposits (or withdrawals). Let k be 
the annual deposit amount, and suppose that an installment is deposited after 
every time interval At. The value of the investment at the time f + At is then 
given by 

S(t + At) = S(t) + (rAt)S(t) + kAt, (2.14) 



where at the end of the time interval At, rAtS(t) is the amount of interest 
credited and kAt is the amount of money deposited (k > 0) or withdrawn 
(k < 0). As a numerical example, if the account held $10,000 at time t, and 
r = 6% per year and k = $12,000 per year, say, and the compounding and 
deposit period is At = 1 month = 1/12 year, then the interest awarded after 
one month is rAtS = (0.06/12) x $10,000 = $50, and the amount deposited is 
kAt = $1000. 



Rearranging the terms of (2.14) to exhibit what will soon become a deriva- 
tive, we have 



S(t + At) - S(t) 
At 



= rS(t) + k. 



The equation for continuous compounding of interest and continuous deposits 
is obtained by taking the limit At — > 0. The resulting differential equation is 



<IS_ 

dt 



rS + k, 



(2.15) 



which can solved with the initial condition <5'(0) = Sq, where Sq is the initial 
capital. We can solve either by separating variables or by using an integrating 
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factor; I solve here by separating variables. Integrating from t = 0 to a final 
time t, 



f s dS 



I S 0 



1 

r 



In 



rS + k 
/ rS + k 
\ rS 0 + k 




rS + k = ( rSo + k)e rt , 
c rSoe rt + ke rt — k 

; 

r 

S = S 0 e rt + ^e rt (1 - e~ rt ) , 



(2.16) 



where the first term on the right-hand side of ( |2.16[ ) comes from the initial 
invested capital, and the second term comes from the deposits (or withdrawals). 
Evidently, compounding results in the exponential growth of an investment. 

As a practical example, we can analyze a simple retirement plan. It is 
easiest to assume that all amounts and returns are in real dollars (adjusted for 
inflation) . Suppose a 25 year-old plans to set aside a fixed amount every year of 
his/her working life, invests at a real return of 6%, and retires at age 65. How 
much must he/she invest each year to have $8,000,000 at retirement? We need 
to solve (2.16) for k using t = 40 years, S(t) = $8,000,000, Sq = 0, and r = 0.06 
per year. We have 



k = 



r S(t) 

e rt - 1 ’ 



0.06 x 8,000,000 

e 0. 06x40 _ l ’ 



= $47,889 year -1 . 



To have saved approximately one million US$ at retirement, the worker would 
need to save about HK$50,000 per year over his/her working life. Note that the 
amount saved over the worker’s life is approximately 40 x $50,000 = $2,000,000, 
while the amount earned on the investment (at the assumed 6% real return) is 
approximately $8,000,000 — $2,000,000 = $6,000,000. The amount earned from 
the investment is about 3 x the amount saved, even with the modest real return 
of 6%. Sound investment planning is well worth the effort. 



2.4.2 Chemical reactions 



Suppose that two chemicals A and B react to form a product C, which we write 
as 

A + B AC, 

where k is called the rate constant of the reaction. For simplicity, we will use 
the same symbol C, say, to refer to both the chemical C and its concentration. 
The law of mass action says that dC/dt is proportional to the product of the 
concentrations A and B, with proportionality constant k ; that is, 



dC 

dt 



= kAB. 



(2.17) 
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Similarly, the law of mass action enables us to write equations for the time- 
derivatives of the reactant concentrations A and B: 



— = - kAB , — = —kAB. 

dt dt 



(2.18) 



The ode given by (2.17) can be solved analytically using conservation laws. We 
assume that An and Bq are the initial concentrations of the reactants, and that 
no product is initially present. From (2.17) and (2.18), 



s<- 1+c > = ° 

5< b + c > = ° 



A + C = A 0 , 
B + C = B 0 . 



Using these conservation laws, (2.17) becomes 

dC 

— = k(A 0 — C)(B 0 — C), (7(0) = 0, 

dt 

which is a nonlinear equation that may be integrated by separating variables. 
Separating and integrating, we obtain 



r C 



dC 



I o (A 0 -C)(B 0 -C) 



= k 



dt 



/o 



= kt. 



(2.19) 



The remaining integral can be done using the method of partial fractions. We 
write 

1 a b 

(A 0 - C)(B 0 - C) = 7L 0 - C + B 0 - C' ( ‘ } 

The cover-up method is the simplest method to determine the unknown coeffi- 
cients a and b. To determine a, we multiply both sides of (2.20) by A 0 — C and 
set C = A 0 to find 

1 

a ~ Bn - An' 



Similarly, to determine 6, we multiply both sides of (2.20) by Bq — C and set 
C = B 0 to find 

b= 1 



Therefore, 



1 



^4o — Bn 



1 



(4) - C)(Bo - C ) B q - A 0 \Ao - C B 0 -C 



1 



1 



and the remaining integral of (2.19) becomes (using C < A 0 ,B 0 ) 



dC 



1 



dC 



dC 



lo ( Ao-CHBo-C ) 



Bq ~ Aq 
1 

Bq — Aq 
1 



— In 



A 0 -C 

A 0 -C 



Bn - C 



In 



An 

Aq(B 0 - C ) 



In 



B 0 -C 

Bn 



Bq — Aq \i?o(4 — C) 
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Using this integral in (2.191, multiplying by ( B 0 — A 0 ) and exponentiating, we 
obtain 



Ap(B 0 C) _ (B 0 -A 0 )kt 
Bo(A 0 -C) 



Solving for C, we finally obtain 



C(t) — AqBq 



e (B 0 -A 0 )kt _ y 
B 0 e( B °~ A O kt - A 0 ’ 



which appears to be a complicated expression, but has the simple limits 



lim C(t) 

t— >00 



Ao if A 0 < B 0 , 
Bo if Bo < Aq, 



= min(A 0 , Bq). 



Hence, the reaction stops after one of the reactants is depleted; and the final 
concentration of product is equal to the initial concentration of the depleted 
reactant. 



2.4.3 The terminal velocity of a falling mass 

view tutorial 

Using Newton’s law, we model a mass m free falling under gravity but with 
air resistance. We assume that the force of air resistance is proportional to the 
speed of the mass and opposes the direction of motion. We define the z-axis to 
point in the upward direction, opposite the force of gravity. Near the surface 
of the Earth, the force of gravity is approximately constant and is given by 
—mg, with g = 9.8 m/s 2 the usual gravitational acceleration. The force of air 
resistance is modeled by —kv, where v is the vertical velocity of the mass and 
k is a positive constant. When the mass is falling, v < 0 and the force of air 
resistance is positive, pointing upward and opposing the motion. The total force 
on the mass is therefore given by F = —mg — kv. With F = ma and a = dv/dt, 
we obtain the differential equation 



dv 

m— = —mq — kv. 
dt y 



(2.21) 



The terminal velocity Uqo of the mass is defined as the asymptotic velocity after 
air resistance balances the gravitational force. When the mass is at terminal 
velocity, dv/dt = 0 so that 



^OO 



mg 

X* 



( 2 . 22 ) 



The approach to the terminal velocity of a mass initially at rest is obtained by 
solving (2.21) with initial condition w(0) = 0. The equation is both linear and 
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separable, and I solve by separating variables: 



TO 



r dv 

' 0 mg + kv 



= - dt, 



to ( mg + kv 

— In ' 



1 + 



V m 9 

kv 



= ~t, 



— g -kt/r 



V = — - 



mg 

mg 



^1 - . 



Therefore, v = Vqo (l — e _fct//m ), and v approaches Voo as the exponential term 
decays to zero. 

As an example, a skydiver of mass to = 100 kg with his parachute closed 
may have a terminal velocity of 200 km/hr. With 



g = (9.8m/s )(10 d km/m)(60 s/min) 2 (60 min/hr) 2 = 127, 008km/hr 



one obtains from (2.22), k = 63, 504 kg/hr. One-half of the terminal velocity 
for free-fall (100 km/hr) is therefore attained when (1 — e -kt/m ) _ jy2, or t = 
77T. In 2/fc ss 4sec. Approximately 95% of the terminal velocity (190 km/hr ) is 
attained after 17 sec. 



2.4.4 Escape velocity 



view tutorial 

An interesting physical problem is to find the smallest initial velocity for a 
mass on the Earth’s surface to escape from the Earth’s gravitational field, the 
so-called escape velocity. Newton’s law of universal gravitation asserts that the 
gravitational force between two massive bodies is proportional to the product of 
the two masses and inversely proportional to the square of the distance between 
them. For a mass m a position x above the surface of the Earth, the force on 
the mass is given by 



F = —G 



Mm 

(. R + x ) 2 ’ 



where M and R are the mass and radius of the Earth and G is the gravitational 
constant. The minus sign means the force on the mass to points in the direction 
of decreasing x. The approximately constant acceleration g on the Earth’s 
surface corresponds to the absolute value of F/m when x = 0: 



9 = 



GM 

~W' 



and g « 9.8 m/s 2 . Newton’s law F = ma for the mass to is thus given by 



d 2 x _ GM 
dt 2 ( R + x) 2 
9 

(1 + x/R) 2 ’ 



(2.23) 



where the radius of the Earth is known to be R ss 6350 km. 
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A useful trick allows us to solve this second-order differential equation as 
a first-order equation. First, note that d 2 x/dt 2 = dv/dt. If we write v(t) = 
v(x(t)) — considering the velocity of the mass m to be a function of its distance 
above the Earth -we have using the chain rule 



dv dv dx 

dt dx dt 

dv 

= V ~dx' 



where we have used v = dx/dt. Therefore, (2.23) becomes the first-order ode 



dv _ g 

V dx (1 +x/R) 2 ' 



which may be solved assuming an initial velocity v(x = 0) = Vq when the mass 
is shot vertically from the Earth’s surface. Separating variables and integrating, 
we obtain 

jC-- jT™- 

The left integral is |(i> 2 — Vq), and the right integral can be performed using 
the substitution u = 1 + x/R, dn = dx/R: 



dx 



(1 + x/R)‘ 



= R 



R 

u 

= R — 
a 

Rx 



r 1+x/R du 

i u 2 

- l+x/R 



R 2 



R 



Therefore, 



*(« 2 -«§) = - 



R 



gRx 



x + R 

which when multiplied by m is an expression of the conservation of energy (the 
change of the kinetic energy of the mass is equal to the change in the potential 
energy). Solving for v 2 , 



v = v n - 



2 gRx 



x + R 

The escape velocity is defined as the minimum initial velocity Vq such that 
the mass can escape to infinity. Therefore, Vq = v escape when v — ► 0 as x — > oo. 
Taking this limit, we have 



2 r 2gRx 

v = lim 

eSCa Pe X + R 



= 2gR. 



With R ps 6350 km and g = 127, 008 km/hr 2 , we determine u eS ca P e = gR ~ 
40 000 km/hr. In comparison, the muzzle velocity of a modern high-performance 
rifle is 4300 km/hr, almost an order of magnitude slower. 
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2.4.5 The logistic equation 

Let N = N(t ) be the size of a population at time t and let r be the growth 
rate. The Malthusian growth model (Thomas Malthus, 1766-1834), similar to 
a compound interest model, is given by 



dN 

— — = rN. 
dt 



Under a Malthusian growth model, a population grows exponentially like 



N(t) = N 0 e rt , 



where No is the initial population size. However, when the population growth 
is constrained by limited resources, a heuristic modification to the Malthusian 
growth model results in the Verhulst equation, 



dN „ A N 

— = rN 1 

dt \ K 



(2.24) 



where K is called the carrying capacity of the environment. Making (2.24) 
dimensionless using r = rt and x = N/K leads to the logistic equation, 



dx 

dr 



x{\ — x), 



where we may assume the initial condition :r(0) = xq > 0. Separating variables 
and integrating 




dx 

x{\ — x) 




dr. 



The integral on the left-hand-side can be done using the method of partial 
fractions: 



1 a b 

x{l — x) x 1 — x 

a + (b — a)x 
x(l — x ) ’ 

and equating the coefficients of the numerators proportional to x° and x 1 , we 
have a = b = 1. Therefore, 



dx 

c(l — x) 



' X 0 



dx 

x 



dx 



xo l 1 - x ) 

= 111 — - In 

x 0 l-Xo 



= In 



x(l - Xq) 
Xq( 1 - x) 



= r. 
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Solving for x, we first exponentiate both sides and then isolate x : 

x(l - x 0 ) _ T 
cr 0 (l - x) 

x{l — Xo) = xoe T — xxoe T , 
x(l — x 0 + x 0 e T ) = x 0 e T , 
x 0 

x = . 

Xo + (1 - xo)e~ T ' 



(2.25) 



We observe that for xq > 0, we have lim T ^. (X) x(t) = 1, corresponding to 



lim N (t) = K. 

t—¥ OO 



The population, therefore, grows in size until it reaches the carrying capacity of 
its environment. 
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Chapter 3 

Second-order linear 
differential equations with 
constant coefficients 



Reference: Boyce and DiPrima, Chapter 3 

The general second-order linear differential equation with independent variable 
t and dependent variable x = x(t) is given by 

x + p(t)x + q(t)x = g(t), (3.1) 

where we have used the standard physics notation x = dx/dt and x = d 2 x/dt 2 . 
A unique solution of © requires initial values x(to) = xo and x(t o) = uq. 
The equation with constant coefficients — on which we will devote considerable 
effort — assumes that p(t) and q(t) are constants, independent of time. The 
second-order linear ode is said to be homogeneous if g(t) = 0. 



3.1 The Euler method 



view tutorial 



In general, (3.1) cannot be solved analytically, and we begin by deriving an 



algorithm for numerical solution. We write the second-order ode as a pair of 
first-order odes by defining u = x, and writing the first-order system as 



where 



x = u, 

u = f{t,x,u ), 

f(t, x, u) = -p(t)u - q(t) x + g(t). 



(3.2) 

(3.3) 



The first ode, (3.2), gives the slope of the tangent line to the curve x = x(t); 
the second ode, (3.3), gives the slope of the tangent line to the curve u = u(t). 



Beginning at the initial values (, x,u ) = (xq,uo) at the time t = to, we move 
along the tangent lines to determine x\ = x(t o + At) and u\ = u(to + At): 



aq = £o + Atzto, 

mi = u 0 + A t/(t 0 , x 0 , m 0 ). 
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The values X\ and zq at the time t\ = t 0 + At are then used as new initial values 
to inarch the solution forward to time t 2 = ti + At. As long as f(t,x,u) is a 
well-behaved function, the numerical solution converges to the unique solution 
of the ode as At — >• 0. 



3.2 The principle of superposition 

view tutorial 

Consider the second-order linear homogeneous ode: 



x + p(t)x + q(t)x = 0; 



(3.4) 



and suppose that x = X-, (t) and x = X l 2 (t) are solutions to (3.4 ) . We consider 
a linear combination of A'i and X 2 by letting 



X(t) = c 1 X 1 (t)+c 2 X 2 (t), (3.5) 



with ci and c 2 constants. The principle of superposition states that x = X{t) 
is also a solution of (|3.4|) . To prove this, we compute 



X + pX + qX = c±Xi + c 2 X 2 + p fciAi + c 2 X 2 ^ + q (ciAi + c 2 X 2 ) 

= ci (Xi + pX i + qXi'j + c 2 (X 2 + pA 2 + qX^ 

= ci x 0 + c 2 x 0 
= 0 , 



since X\ and AT 2 were assumed to be solutions of (3.4). We have therefore shown 
that any linear combination of solutions to the second-order linear homogeneous 
ode is also a solution. 



3.3 The Wronskian 



view tutorial 



Suppose that having determined that two solutions of (3.4) are x = A'i (t) and 



x = X 2 (t), we attempt to write the general solution to (3.4) as (3.5). We must 



then ask whether this general solution will be able to satisfy the two initial 
conditions given by 

x{t 0 ) = x 0 , x(t 0 ) = uq . (3.6) 



Applying these initial conditions to (3.5), we obtain 



CiX^tg) + c 2 AT 2 (t 0 ) — Xq, 
Cl Xi (to) + c 2 A 2 (to) = Uq, 



(3.7) 



which is observed to be a system of two linear equations for the two unknowns 



Ci and c 2 . Solution of (3.7) by standard methods results in 



Ci = 



^oA 2 (to) ~ A 2 (to) 

W ’ 



c 2 = 



wo A’i (to) — ^oAi(to) 

W" 
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where W is called the Wronskian and is given by 

W = X^toJMto) - X 1 (to)X 2 (to). (3.8) 

Evidently, the Wronskian must not be equal to zero ( W ^ 0) for a solution to 
exist. 

For examples, the two solutions 

Xi{t) = sinwt, X 2 {t) = Asinwi, 

have a zero Wronskian at t = to, as can be shown by computing 

W = (sinwto) (Aucoswfo) — (ujcosuto) (Asinwto) 

= 0 ; 

while the two solutions 



X\ (t) = sin cut, X 2 (t) = cos wt, 
with u/0, have a nonzero Wronskian at t = to, 

W = (sinwfo) (— wsinwto) — (wcoswto) (coswfo) 

= —uj. 

When the Wronskian is not equal to zero, we say that the two solutions 
X\ (t) and X 2 (t) are linearly independent. The concept of linear independence 
is borrowed from linear algebra, and indeed, the set of all functions that satisfy 
(|3.4|) can be shown to form a two-dimensional vector space. 



3.4 Second-order linear homogeneous ode with 
constant coefficients 



view tutorial 

We now study solutions of the homogeneous, constant coefficient ode, written 
as 

ax + bx + cx = 0, (3-9) 

with a, 6, and c constants. Such an equation arises for the charge on a capacitor 
in an unpowered RLC electrical circuit, or for the position of a freely-oscillating 
frictional mass on a spring. Our solution method finds two linearly independent 



solutions to (3.9), multiplies each of these solutions by a constant, and adds 



them. The two free constants can then be used to satisfy whatever initial 
conditions are given. 

Because of the differential properties of the exponential function, a natural 
ansatz, or educated guess, for the form of the solution to (3.9) is x = e rt , where 



r is a constant to be determined. Successive differentiation results in x = re rt 



and x = r 2 e rt , and substitution into (3.9) yields 



ar 2 e rt + bre rt + ce rt = 0. 



(3.10) 



Our choice of exponential function is now rewarded by the explicit cancelation 



of e rt in (3.10). The result is a quadratic equation for the unknown constant r: 



ar 2 + br + c = 0. 



(3.11) 
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Our ansatz has thus converted a differential equation into an algebraic equation. 



Equation ( 3.11 ) is called the characteristic equation of ( 3.9 ). Using the quadratic 



formula, the two solutions of the characteristic equation (3.11 1 are given by 



r± = 



2a 



b± \Jb 2 — 4ac^ . 



There are three cases to consider: (1) if b 2 — 4ac > 0, then the two roots are 
distinct and real; (2) if b 2 — lac < 0, then the two roots are distinct and complex 
conjugates of each other; (3) if b 2 — 4 ac = 0, then the two roots are degenerate 
and there is only one real root. We will consider these three cases in turn. 



3.4.1 Real, distinct roots 



When r_|_ ^ r_ are real roots, then the general solution to (3.9) can be written 
as a linear superposition of the two solutions e r+t and e r_t ; that is, 

x(t) = Cie r+t + C2e r ~ t . 



The unknown constants c\ and C 2 can then be determined by the initial condi- 
tions a; (to) = Xq and a; (to) = uo- We now present two examples. 



Example 1: Solve x + 5± + 6a: = 0 with x(O) = 2, i(O) = 3, and find the 
maximum value attained by x. 

view tutorial 

We take as our ansatz x = e rt and obtain the characteristic equation 

r 2 + 5r + 6 = 0, 



which factors to 

(r + 3)(r + 2) = 0. 

The general solution to the ode is thus 

a ;(t) = ciW 2t + C2e~ 3t . 

The solution for x obtained by differentiation is 

x(t) — — 2cie -2t — 3c2e _3t . 

Use of the initial conditions then results in two equations for the two unknown 
constant C\ and C 2 : 



Ci + C 2 — 2, 

— 2ci — 3c2 = 3. 

Adding three times the first equation to the second equation yields C\ = 9; and 
the first equation then yields C 2 = 2 — ci = — 7. Therefore, the unique solution 
that satisfies both the ode and the initial conditions is 

a :(t) = 9e~ 2t - 7e~ 3t 

= 9 '- 2, ( '-H- 
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Note that although both exponential terms decay in time, their sum increases 
initially since x(0) > 0. The maximum value of x occurs at the time t m when 
x = 0, or 

t m = In (7/6) . 

The maximum x m = x(t m ) is then determined to be 

x m = 108/49. 



Example 2: Solve x — x = 0 with x(0) = Xo, x(0) = uq- 

Again our ansatz is x = e rt , and we obtain the characteristic equation 

r 2 - 1 = 0, 

with solution r± = ±1. Therefore, the general solution for x is 

x(t) = cie* + c 2 e~ t , 

and the derivative satisfies 

x(t) = cie* — c 2 e~ l . 

Initial conditions are satisfied when 



Cl + c 2 — Xo, 

Cl — c 2 = Uo- 

Adding and subtracting these equations, we determine 

ci = ^ {xq +u 0 ), c 2 = ^ (x 0 - w 0 ) , 
so that after rearranging terms 




The terms in parentheses are the usual definitions of the hyperbolic cosine and 
sine functions; that is, 

. e t + e _t . , e* — e _t 

cosh t = , smh t = . 

2 ’ 2 

Our solution can therefore be rewritten as 

x(t) = Xq cosh t + uq sinh t. 



Note that the relationships between the trigonometric functions and the complex 
exponentials were given by 



cos t = 




sint = 




so that 



coshi = cos it, sinht = —i sin it. 



Also note that the hyperbolic trigonometric functions satisfy the differential 
equations 

d . d 

— smh t = cosh t, — cosh t = smh t, 
dt dt 

which though similar to the differential equations satisfied by the more com- 
monly used trigonometric functions, is absent a minus sign. 
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3.4.2 Complex conjugate, distinct roots 

view tuto rial 

We now consider a characteristic equation (3.11) with b 2 —4 ac < 0, so the roots 
occur as complex conjugate pairs. With 



fi = ^-\J4 ac-b 2 , 
2 a 2 a 



the two roots of the characteristic equation are A + ifi and A — ifi. We have 
thus found the following two complex exponential solutions to the differential 
equation: 

Z 1 = e xt e i » t , Z 2 = e xt e~ iflt . 

Applying the principle of superposition, any linear combination of Z\ and Z 2 
is also a solution to the second-order ode, and we can form two different linear 
combinations that are real. They are given by 



and 



Xi (t) 



Z\ + z 2 
2 




e xt cos fit, 



x 2 (t) 



Z i — z 2 
2 i 




e xt sin fit.. 



Having found the two real solutions Xi(f) and A 2 (f), we can then apply the 
principle of superposition again to determine the general solution x(t) to be 

x(t) = e xt (A cos fit + B sin fit) . (3-12) 



It is best to memorize this result. The real part of the roots of the characteristic 
equation goes into the exponential function; the imaginary part goes into the 
argument of cosine and sine. 



Example 1: Solve x + x = 0 with x(0) = Xo and x(0) = ito 

view tutorial 

The characteristic equation is 

r 2 + 1 = 0, 



with roots 



r± = ±i. 



The general solution of the ode is therefore 



x(t) = A cost + B sint. 
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The derivative is 

x(t) = — A sin t + B cos t. 

Applying the initial conditions: 

x(0) = A = Xo, x(0) = B = uq; 
so that the final solution is 

x(t) = Xo cos t + uq sint. 

Recall that we wrote the analogous solution to the ode x — x = 0 as x(t) = 
xq cosh t + Uq sinlr t. 



Example 2: Solve x + x + x = 0 with a:(0) = 1 and i(0) = 0. 
The characteristic equation is 

r 2 + r + 1 = 0, 

with roots 

1 , V3 

r± = — - ± i —— . 

2 2 

The general solution of the ode is therefore 



x(t) = e 2 4 [ A cos + -Bsin ] . 



The derivative is 

1 



^ ( Vs . Vs 

x(t) = — -e 2 Acos — - 1 + Bsm — t 



\/S ( A . VS „ \ 

H — — e 2 * —Asm — — t + B cos — — t . 



Applying the initial conditions x(0) = 1 and i:(0) = 0: 

A = 1, 

_A A+ V^ b = 0; 



or 



Therefore, 






. , _i t ( Vs Vs . Vs \ 

x(t) = e 2 * cos —t + — sin —t . 
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3.4.3 Repeated roots 

view tutorial 

Finally, we consider the characteristic equation, 

ar 2 + br + c = 0, 

with b 2 — 4 ac = 0. The degenerate root is then given by 

b 

r ~ _ 2a’ 

yielding only a single solution to the ode: 

xi{t) = exp ( 3 - 13 ) 

To satisfy two initial conditions, a second independent solution must be found, 
and apparently this second solution is not of the form of our ansatz x = exp (rt) . 
One method to determine this missing second solution is to try the ansatz 



x(t) = y(t)xi(t), 



(3.14) 



where y(t) is an unknown function that satisfies the differential equation ob- 
tained by substituting (3.14| into (3.9 1. This general technique is called the 



reduction of order method and enables one to find a second solution of a ho- 
mogeneous linear differential equation if one solution is known. If the original 
differential equation is of order n, the differential equation for y = y(t ) reduces 
to one of order n — 1 . 

Here, however, we will determine this missing second solution through a 
limiting process. We start with the solution obtained for complex roots of the 
characteristic equation, and then arrive at the solution obtained for degenerate 
roots by taking the limit y — >■ 0. 



Now, the general solution for complex roots was given by (3.12), and to 
properly limit this solution as y — > 0 requires first satisfying the specific initial 
conditions x(0) = xq and x(0) = uq. Solving for A and H, the general solution 



given by (3.12) becomes the specific solution 



, . \t ( u 0 — Axo 

x(t; y) = e x 0 cos/if-| sin yt 

V M 

Here, we have written x = x(t; y) to show explicitly that x depends on y. 
Taking the limit as y — > 0, and using lim^o y -1 sin yt = t , we have 

lim x(t; y) = e xt (xq + (uq — Xxq ) t) . 

f-L — ^0 V 7 

The second solution is observed to be a constant, up — Axo, times t times the 
first solution, e At . Our general solution to the ode (3.9) when b 2 — 4ac = 0 can 
therefore be written in the form 



x{t) = (ci + c 2 f)e rt , 

where r is the repeated root of the characteristic equation. The main result to 
be remembered is that for the case of repeated roots, the second solution is t 
times the first solution. 
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Example: Solve x + 2x + x = 0 with a;(0) = 1 and i(0) = 0. 

The characteristic equation is 

r 2 + 2r + 1 = (r + l) 2 
= 0 , 

which has a repeated root given by r = — 1. Therefore, the general solution to 
the ode is 

x(t) = Cie _t + C 2 te~ t , 

with derivative 

x(t) = —Cie~ t + c 2 e~ l — c 2 te~ t . 

Applying the initial conditions: 

Ci = 1, 

-ci + c 2 = 0, 

so that Ci = c 2 = 1. Therefore, the solution is 

x(t) = (1 + i)e _t . 



3.5 Second-order linear inhomogeneous ode 



We now consider the general second-order linear inhomogeneous ode (3.1 1 : 

x +p(t)x + q(t)x = g(t), (3.15) 

with initial conditions x(to ) = Xq and x(to) = uo. There is a three-step solution 
method when the inhomogeneous term g(t ) 7^ 0. (i) Find the general solution 
of the homogeneous equation 



x + p(t)x + q{t)x = 0. 

Let us denote the homogeneous solution by 

Xh(t) = CiXi(t) + c 2 X 2 (t), 



(3.16) 



where Xi and A” 2 are linearly independent solutions of (3.161, and Ci and c 2 



are as yet undetermined constants, (ii) Find any particular solution x p of the 



inhomogeneous equation (3.15). A particular solution is readily found when 



p(t) and q(t) are constants, and when g(t) is a combination of polynomials, 



exponentials, sines and cosines, (iii) Write the general solution of (3.151 as the 
sum of the homogeneous and particular solutions, 



x(t) = x h (t) +x p (t), 



(3.17) 



and apply the initial conditions to determine the constants c\ and c 2 . Note that 



because of the linearity of (3.15) 



x+px + qx = -i^{ x h + x p ) +p—(x h + x p ) + q(x h + x p ) 

= ( Xh +pXh + qxh) + (■ Xp +px p + qxp) 

= 0 + g 
= 9, 



